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We report the results of a systematic investigation of the phase diagram of the iron-based su- 
perconductor, Bai_ a ,K a! Fe2As2, from x — to x — 1.0 using high resolution neutron and x-ray 
diffraction and magnetization measurements. The polycrystalline samples were prepared with an 
estimated compositional variation of Ax < 0.01, allowing a more precise estimate of the phase 
boundaries than reported so far. At room temperature, Bai_ a; K I Fe2As2 crystallizes in a tetragonal 
structure with the space group symmetry of I4/mmm, but at low doping, the samples undergo 
a coincident first-order structural and magnetic phase transition to an orthorhombic (O) struc- 
ture with space group Fmmm and a striped antiferromagnet (AF) with space group F c mm'm'. 
The transition temperature falls from a maximum of 139 K in the undoped compound to K at 
x — 0.252, with a critical exponent as a function of doping of 0.25(2) and 0.12(1) for the structural 
and magnetic order parameters, respectively. The onset of superconductivity occurs at a critical 
concentration of x = 0.130(3) and the superconducting transition temperature grows linearly with 
x until it crosses the AF/O phase boundary. Below this concentration, there is microscopic phase 
coexistence of the AF/O and superconducting order parameters, although a slight suppression of 
the AF/O order is evidence that the phases are competing. At higher doping, superconductivity has 
a maximum T c of 38 K at x = 0.4 falling to 3 K at x = 1.0. We discuss reasons for the suppression 
of the spin-density-wave order and the electron-hole asymmetry in the phase diagram. 



I. INTRODUCTION 

There is now an extensive body of research into the ori- 
gin of superconductivity in the iron-based superconduc- 
tors demonstrating the importance of the subtle interplay 
of their electronic properties with crystalline structure 
[1,2]. These compounds all contain a common structural 
motif, namely a square planar net of iron atoms tetra- 
hcdrally coordinated with pnictogens or chalcogens pro- 
ducing Fe 2 A 2 layers, in which X — As or Se/Te in the 
highest T c compounds. They are separated by buffer lay- 
ers comprising, for example, rare earth oxides in the so- 
called 1111 systems, such as those based on LaFeAsO, or 
alkaline earths in the so-called 122 systems, such as those 
based on BaFe 2 As 2 . Like the cuprate superconductors, 
the buffer layers can act as charge reservoirs, controlling 
the carrier concentration and inducing superconductiv- 
ity in the iron planes by the introduction of aliovalent 
dopants, e.g., LaFeAsOi-zF^ [3,4] and Bai_ a; K :E Fe 2 As 2 
[5,6], but it is also possible to dope the pnictogen or 
chalcogen sites, e.g., BaFe 2 As 2 _ 2: P 2: [7,8], or the iron 
planes themselves by substituting other transition metal 
ions, e.g., BaFe 2 _ x Co x As2 and BaFe 2 _ a; Ni :E As2 [9,10]. 

In this article, we report on a systematic investigation 
of the phase diagram of Bai_ :E K x Fe 2 As2 from x = to 
x = 1.0 using high resolution neutron and x-ray diffrac- 
tion combined with bulk characterization. In spite of the 
diversity of doping strategies employed to modify the su- 
perconducting properties of the iron-based superconduc- 
tors, their phase diagrams show remarkable similarities. 
There is typically an undoped "parent" compound that is 



antifcrromagnetic rather than superconducting[l]. These 
are fully compensated metals, whose electronic properties 
are dominated by multiple iron-derived cZ-bands near the 
Fermi level with approximately equal concentrations of 
hole and electron carriers [11]. The Fermi surfaces con- 
sist of quasi-two-dimensional cylinders, with two or three 
hole pockets at the Brillouin zone centers and two elec- 
tron pockets at the M-points on the zone boundaries, 
i. e. , along the direction of the nearest-neighbor iron-iron 
bonds [12]. All the Fermi surface pockets have similar 
radii in the undoped compounds, making their electronic 
structure particularly susceptible to magnetic instabili- 
ties resulting from a nesting of the disconnected hole and 
electron Fermi surfaces [13]. Since the nesting wavevec- 
tor corresponds to the antifcrromagnetic wavevector ob- 
served by neutron diffraction [14-16], it is plausible that 
the magnetism can be explained by a purely itinerant 
model of spin density waves, and ab initio density func- 
tional theory does indeed predict the correct magnetic 
structure [17], although there is an ongoing debate about 
the strength of electron correlations [18,19]. 

The magnetic structure breaks the tetragonal sym- 
metry with an in-plane wavevector of cither (0,7r) or 
(7r,0) in the unfolded Brillouin zone with one iron atom 
per unit cell. With a finite magnetoelastic coupling, 
this would induce an orthorhombic structural transition, 
which is usually observed to occur at the same tem- 
perature as magnetic order in the parent compounds 
[20,21]. However, the addition of both hole and electron 
charge carriers through chemical substitution suppresses 
both transitions. A variety of scenarios are possible in 
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FIG. 1: (Color online) Structure of BaFe2As2, which crystal- 
lizes in a tetragonal ThCr2Si2-type structure with the space 
group symmetry of 14/mmm. Potassium substitutes onto the 
barium sites. 

Ginzburg-Landau treatments of the magnetoelastic cou- 
pling [22,23]. The two phase transitions could be first- 
or second-order and occur simultaneously or separately. 
In most of the iron-based compounds, the two transi- 
tion temperatures split with doping [24] and there is a 
report of a split transition in Bai_ a; K x Fe2As2 as well 
[25]. However, this is inconsistent with our previously 
reported neutron diffraction data, which shows unam- 
biguously that they are coincident and first-order for all 
x before they are both suppressed at x < 0.3 [6]. Unusu- 
ally, the two order parameters, magnetic and structural, 
are proportional to each other, apparently indicating a 
biquadratic coupling that is usually only observed at a 
tetracritical point [26,27], not over an extended range of 
compositions. Possible explanations for this observation 
will be discussed in the conclusions. 

Superconductivity emerges before the complete sup- 
pression of the antiferromagnetic/orthorhombic (AF/O) 
phase and coexists at low-doping levels. The nature of 
the competition between AF/O order and superconduc- 
tivity is a central question in understanding iron-based 
superconductivity [28,29]. There were earlier reports 
based on local probes that, in Bai_ a; K. c Fe2As2, the coex- 
istence region is characterized by a mesoscopic phase sep- 
aration into AF/O and superconducting droplets [30,31], 
but our previously reported diffraction data are only con- 
sistent with a microscopic phase coexistence [6], a con- 
clusion since supported by muon spin rotation (/iSR) ex- 
periments [32], suggesting that the earlier reports may 
be due to compositional fluctuations within the samples. 

One of the main reasons for studying Bai_ ;E K. E Fe2As2 
is that superconductivity extends up to much higher 
hole-doping levels, with 0.5 holes/Fe atom, than in the 
electron-doped superconductors produced by transition 
metal substitutions. In the case of BaFe2-^Co x As2, su- 



perconductivity vanishes at only 0.12 electrons/Fe atom 
[10]. Furthermore, the maximum T c with hole-doping 
is 38 K, significantly higher than the maximum T c of 
~ 25 K obtained with electron doping. This electron- 
hole asymmetry in the phase diagram has been attributed 
to an enhanced Fermi surface nesting in the hole-doped 
compounds, consistent with angle-resolved photoemis- 
sion spectroscopy (ARPES) data and band structure cal- 
culations [33] . This explanation is also supported by the 
evolution of resonant spin excitations, which become in- 
commensurate due to the mismatch in hole and electron 
Fermi surface volumes when T c starts to fall [34] . On the 
other hand, there is also a strong correlation between T c 
and internal structural parameters such as the Fe-As-Fe 
bond angles [35,36]. These are known to have an influ- 
ence on the band structure and the degree of moment 
localization, but their role in optimizing superconductiv- 
ity and the implications for the gap symmetry is a matter 
of debate [37,38]. 

There have been two previous reports of the doping de- 
pendence of this series in addition to our own brief report, 
which are all in qualitative agreement [6,36,39]. At room 
temperature, all members of the Bai_ :c K 2 ;Fe2As2 series 
crystallize in a tetragonal structure with the space group 
symmetry of IA/mmm (Fig. 1), while low-doped samples 
also exhibit a low temperature phase transition to an or- 
thorhombic structure with space group Fmmm [5]. Su- 
perconducting samples at higher doping have a maximum 
T c of 38 K and remain tetragonal at all measured temper- 
atures down to 1.7 K. However, there are significant dis- 
crepancies in the published reports concerning the criti- 
cal dopant concentrations defining the onset of supercon- 
ductivity and the suppression of the AF/O phase, with 
the latter varying from x ~ 0.3 [36] to ~ 0.4 [39]. As 
already mentioned, there have also been disagreements 
about the nature of the competition between the three 
ordered phases at low doping. We believe that these dis- 
crepancies are due to uncertainties in the actual compo- 
sition of the synthesized samples, since it is well-known 
that potassium is particularly volatile. Controlling the 
inhomogeneity to within acceptable limits in order to im- 
prove the accuracy of the various phase boundaries has 
been a key goal of this work and we estimate that we 
have been able to make samples in which Ax < 0.01. We 
have performed neutron and x-ray diffraction studies of 
the magnetic and structural order using high-resolution 
powder diffractometers so that the systematic variation 
of the lattice parameters and internal structural param- 
eters can be used to estimate the degree of uncertainty 
in the average composition and its variation within the 
samples. 

In this article, we present the results of Rietveld refine- 
ments for the entire series and use this analysis along with 
bulk measurements to produce a comprehensive magnetic 
and structural phase diagram that provides insight into 
the nature of the phase competition that underlies iron- 
based superconductivity. Our results show that there 
is a steeper decrease in T c and hence a narrower region 
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TABLE I: Structural and magnetic phase diagram of Bai- 2 ,K a: Fe2 As2. The nominal value of x represents the starting compo- 
sition given by the Ba/Fe ratio. The fitted value is determined by smoothing the variation in the a-axis lattice parameter from 
0-1 < x < 0.25 using a power law function. The first column is the value of x used in the text and figures. The superconducting 
transition temperature is determined from magnetization measurements. Tm and T 3 are determined by magnetization and neu- 
tron diffraction data as described in the text. The magnetic moments and orthorhombic order parameters, S — (a — b)/(a + b), 
are determined from the low-temperature Rietveld refinements. 



of phase coexistence of the AF/O order with supercon- 
ductivity than previous reports. After a description of 
our experimental results, we combine our findings with 
results reported in the literature on the electron-doped 
BaFe2_ a; Co :I ;As2 series in order to elucidate the origin of 
the electron-hole asymmetry in the phase diagram. 

II. EXPERIMENTAL DETAILS 

The synthesis of homogeneous single phase 
Bai_ a; K a ;Fe2As2 samples is known to be particu- 
larly delicate due to unfavorable kinetics, high vapor 
pressures, and a significant difference in the chemical 
reactivity of K and Ba with FeAs that may result in 
stabilizing other binary by-products. For this work, 
the synthesis and properties of our samples were opti- 
mized by the systematic examination of all reasonable 
combinations of reaction parameters, e.g., purity of the 
starting materials, reaction containers, temperature, 
and duration of heating, etc. Our final samples were 
produced according to the following procedure: Handling 
of all materials was performed in a nitrogen-filled glove 
box. Raw materials (BaAs/KAs/Fe2As) were prepared 
by heating elemental mixtures at 400° C, 600° C, and 
850° C, respectively. The stoichiometric mixture of 
these starting precursors for a desired composition was 
thoroughly ground to ensure uniform and homogeneous 



and subsequently annealed at 1050° C in a Nb tube 
sealed in a quartz tube. The closed metal tubes are 
needed to eliminate any chemical loss that may oth- 
erwise result from the evaporation of K and As. A 
large number of high quality samples were synthesized 
covering the full phase diagram < x < 1 with special 
emphasis given to the 0.1 < x < 0.25 range in which 
the potassium content was incremented in very small 
amounts with Ax = 0.025. The deliberate synthesis of 
samples with finely tuned K content was necessary in 
order to carefully investigate the rapid suppression of 
magnetism with increasing K and to elucidate the nature 
of phase coexistence with superconductivity within the 
same sample. Samples with coarse K increments would 
otherwise lead to inconclusive results. 

Initial characterization of the samples was performed 
by x-ray diffraction, magnetization measurements, and 
Inductively Coupled Plasma (ICP) elemental analysis. 
For select samples, neutron powder diffraction experi- 
ments were performed on the HRPD (x = 0, 0.1, 0.125, 
0.15, 0.175, 0.2, 0.21, 0.24, 0.3, 0.5 and 1) and Wish 
diffractometers (x = 0.22 and 0.25) at the ISIS Pulsed 
Neutron and Muon Source (Rutherford Appleton Labo- 
ratory) and on beamline 11-BM, (x = 0.28) at the Ad- 
vanced Photon Source (Argonne National Laboratory). 
The resolution, Ad/d, at 2 Ais 0.001 for HRPD and 0.002 
for Wish. For diffraction experiments, the samples were 
sealed under vacuum for shipment and re-opened just be- 
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fore the measurements in a helium environment to pre- 
vent air exposure. The nuclear and magnetic structures 
together with interatomic bond-lengths and bond-angles 
were determined by the Rietveld refinement technique 
using the comprehensive General Structure Analysis Sys- 
tem (GSAS) software suite [40] and the associated graph- 
ical user interface (EXPGUI) [41]. Minute traces of no 
more than 0.5-3% by weight of FeAs and Fe2 As impurity 
phases were observed in some of the samples, too small 
to affect the analysis. 

Our neutron diffraction results show that the synthe- 
sis methods has reduced the compositional uncertainty 
significantly, which we estimate to be Aa; < 0.01. In 
previously reported phase diagrams, it is not totally 
clear whether the potassium contents were nominal or 
actual measured values. Because of its volatility, it is 
always necessary to add excess potassium, so the even- 
tual stoichiometry is largely governed by the Ba/Fe ra- 
tio. Thus, any perceived discrepancy with other work, 
such as Ref. [39] is probably due to differences in the 
handling and control of the volatile potassium and ar- 
senic constituents. In this article, we have used a num- 
ber of methods to characterize the sample compositions, 
including direct measurements of the stoichiometry from 
ICP analysis. To produce the final compositions used in 
our phase diagrams, we started with the nominal x val- 
ues determined from the starting Ba/Fe ratio and then 
smoothed the variation in the a-axis lattice parameter 
from 0.1 < x < 0.25 using a power law function (see Ta- 
ble 1). In most cases, the agreement with the nominal 
value was better than 0.005, with just two samples re- 
quiring a significant shift of 0.02 and 0.03, respectively. 
In all cases, the adjustments also improved the consis- 
tency of other measurements, such as the variations in 
transition temperatures and order parameters. The first 
column of Table 1 shows the value of x that we have used 
in the text and figure labels, derived by rounding the fit- 
ted values to the nearest 0.005, although the fitted values 
were used in the plots and in numerical analyses of the 
doping dependence, e.g., the fit of T c vs x. For x > 0.28, 
where the precise composition is not so critical, we have 
used the nominal values. 

We are also able to monitor fluctuations in compo- 
sition within a single sample because HRPD has suffi- 
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FIG. 2: (Color online) Variation of lattice constants, a, b and 
c, with temperature for x — 0, 0.1, 0.21 and 0.3. The lattice 
constants, a and b, in the orthorhombic phase are divided by 
V2. 



ciently high resolution that it is sensitive to distributions 
of the lattice parameter and internal strains caused by 
compositional gradients and even particle size broaden- 
ing [42]. The diffraction peaks from the (220) reflec- 
tion shown in Fig. 1 of Ref. [6] show that there is no 
change in the linewidths and lineshapes of the diffraction 
peaks from the undoped compound up to x = 0.24, con- 
sistent with a high-degree of compositional homogeneity 
(Ax < 0.01). The regular spacing between the peak po- 
sitions is in agreement with the fixed steps in x. 



III. RESULTS 
A. Structural Phase Diagram 

In the undoped BaFe2As2 compound, there is a struc- 
tural phase transition at 139 K from the tetragonal 
ThCr2Si2-type structure with the space group symmetry 
I4/mmm to an orthorhombic /3-SrRIi2As2-type structure 
with space group Fmmm [5] . The structure of both the 
tetragonal and orthorhombic phases of Bai_ a; K 2 .Fe2As2 
can be described as a stack of edge-sharing Fe2As2 layers 
separated by layers of (Ba,K) ions (Fig. 1). The (Ba,K) 
ions occupy crystallographic positions that are tetrahe- 
drally coordinated with four arsenic anions. 

With potassium doping, the tetragonal to orthorhom- 
bic structural transition temperature, T s , decreases until 
it is fully suppressed for x > 0.25. In Fig. 2, we show 
the lattice parameters as a function of temperature for 
x = 0, 0.1, 0.21 and 0.3. Below x = 0.3, a significant 
orthorhombic splitting of the basal plane a and b lattice 
parameters is observed. The evolution of the orthorhom- 
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FIG. 3: (Color online) Temperature dependence of the (110) 
peak in the vicinity of structural transition temperature T s . 
The bold curve shows the peak at the estimated T 3 . For 
< x < 0.24, the other curves show the peak at intervals 
of 2K around T s . For x = 0.3, the peak is shown in 20 K 
intervals between 1.7 K and 120 K. 
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FIG. 4: (Color online) Variation of lattice constants and vol- 
ume in Bai_i,K :E Fe2As2 with x at 1.7 K. Solid lines are guide 
to the eye. The lattice constants (a and b) and the volume in 
the orthorhombic phase are divided by V2 and 2, respectively, 
from those for the Fmmm space group. 

bic order parameter defined by 5 — (a — b)/(a + 6) is 
discussed later. Although the transitions appear to be 
continuous, we observed small but sharp volume anoma- 
lies at all the structural phase transitions (Fig. 3 in Ref. 
[6] ) showing that they are weakly first-order in character 
over the entire phase diagram. 

The structural transition temperatures in Table 1 were 
determined by fitting a power law, 5 cx (T S ~T)P /T s , close 
to the transition temperature, yielding exponents mostly 
in the range /3 ~ 0.13 to 0.2. Apart from x = 0, where 
[3 = 0.129(3), in reasonable agreement with Ref. [27], the 
exponents will be modified by compositional fluctuations 
and so are not reliable estimates of the critical behavior. 
The exponents at x = 0.21 and 0.24 were anomalously 
high (/3 = 0.25 and 0.30, respectively), which could re- 
flect a slightly greater degree of compositional variation 
within those samples. As a check on these values of T s , we 
used the peak profiles close to the transition temperature. 
Fig. 3 shows the temperature dependence of the tetrag- 
onal (110) Bragg peak for x = 0, 0.1, 0.15, 0.2, 0.24 and 
0.3 in the vicinity of the structural transition. This reflec- 
tion splits into two (022) and (202) orthorhombic peaks 
below T s . Close to the transition, the two peaks can- 
not be resolved but we can determine T s from the tem- 
perature dependence of the full width at half maximum 
(FWHM) as they merge. In the high-temperature phase, 
the FWHM of ~ 0.0037(3)Ais independent of the com- 
position. These two methods of determining T s agreed 
within the errors. 

Fig. 4 shows the evolution of the lattice parameters 
and unit cell volume at 1.7 K as a function of x. Barium 
substitution by the smaller potassium cations reduces the 
in-plane a- and 6-lattice parameters while significantly 
lengthening the out-of-plane c-axis. However, the c-axis 
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FIG. 5: (Color online) Variation of Fe-Fe and Fe-As bond 
lengths with x at 1.7 K and with temperature for different 
K substitutions. Blue triangles represent the Fe-As bonds. 
Black square and red circle symbols represent the Fe-Fe bond 
lengths merging at T s . Solid lines are guides to the eye. 



enhancement is not large enough to fully compensate for 
the shrinking basal plane axes and the unit cell volume 
gradually decreases in magnitude upon increasing the K 
content until x ~ 0.5. The reason for the non-monotonic 
behavior of the unit cell volume at high dopant levels is 
not understood and will require further investigation. 

The behavior of the Fe-Fe and Fe-As interatomic dis- 
tances (at 1.7 K) are presented in Fig. 5. The Fe-Fe dis- 
tances mimic the in-plane lattice parameters both as a 
function of K content and of temperature. The similarity 
in behavior is explained by the fact that Fe atoms occupy 
special rigid positions along the long edges of the lattice. 
Six As-Fe-As bond-angles can be identified in each FeAs4 
tetrahedron. In the tetragonal I4/mmm structure, these 
angles can be grouped into two independent angles: two 
equivalent angles, ai, and four equivalent smaller ones, 
«2- In the orthorhombic Fmmm structure below Tjv, 
the angle ct\ remains unaffected but the angle a 2 splits 
into two pairs of equivalent angles a' 2 and a 2 with an 
angular separation of ~ 0.5 — 0.6°. A sketch showing the 
different angles is displayed in Fig. 6. As shown in the 
same panel, the angle ol\ increases linearly and contin- 
uously throughout the whole phase diagram. However, 
starting from BaFe2As2, the angles a' 2 and a' 2 ' increase 
with increasing K until reaching a critical composition 
below x ~ 0.3 beyond which the structural transitions 
are suppressed and the two angles merge into a 2 , which 
continues to increase with higher K contents. The refined 
values of the lattice parameters, bond lengths and bond 
angles at 1.7 K are shown in Table 2. 
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5.5636(4) 


5.5476(4) 


13.1027(10) 


404.157(90) 


0.35398(6) 


2.3900(5) 


2.7810(2) 


2.7732(2) 


108.846(17) 


109.073(16) 


110.50(3) 


0.3 


3.9165(2) 
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201.554(20) 


0.35376(4) 


2.3859(3) 


2.7501(1) 




109.615(11) 




109.185(21) 


1 


3.8251(2) 




13.7846(5) 


201.691(40) 


0.35314(4) 


2.3833(4) 


2.7047(1) 




110.855(12) 




106.708(23) 



TABLE II: Results of Rietveld refinements for Bai_ a; K I Fe2As2 from neutron powder diffraction data collected on HRPD at 
1.7K. For x < 0.3, the space group is Fmmm, in which a 7^ b, there are two inequivalent Fe-Fe bond distances and three 
inequivalent As-Fe-As bond angles. For x > 0.3, the space group is 74/mmm, in which a — b, there is one Fe-Fe bond distance 
and two inequivalent As-Fe-As bond angles 

B. Magnetic Phase Diagram 

Neutron powder diffraction data reveal the presence of 
weak magnetic Bragg reflections that appear below the 
structural phase transition for all of the orthorhombic 
samples. The magnetic peaks shown in Fig. 7, located at 
2.45 Aand 3.43 A, were indexed as 121 and 103 in agree- 
ment with the widely reported antiferromagnetic spin 
density wave (SDW) ground state [15,16,39,43]. As with 
the structural transitions, the antiferromagnetic transi- 
tion (Neel) temperatures as a function of doping were 
determined by power-law fits to temperature variation 
of the magnetic moment (see Table 1). These coincide 
with the orthorhombic transition for all values of x. It 
is important to realize that the structural and magnetic 
orders are identified in the same measurements, the first 
from the splitting of the nuclear Bragg peaks and the sec- 
ond by the intensity of the magnetic Bragg peaks. This 
means that the conclusion that the two transitions are 
coincident does not depend on the accuracy of the ther- 
mometry. As shown in Ref. [6] , the two order parameters 
determined in this way are directly proportional to each 
other at all temperatures over the entire phase diagram, 
an unusual result that we will discuss in more detail in 
the discussion. 

As a further check on our data, the values of Tjy de- 
termined by neutron diffraction are in excellent agree- 
ment with those determined by peaks in the tempera- 
ture derivative of the magnetization (Fig. 8 and Table 
1). These magnetization peaks decrease in magnitude be- 
cause of the progressive attenuation of the magnetic sig- 
nal due to increasing K content until they arc no longer 
detected for x > 0.21. 

A full analysis of the magnetic structure was per- 
formed using the allowed subgroup magnetic symmetries 
of Fmmm. All possible models were tested but only the 
magnetic space group F c mm'm' resulted in a proper fit 
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FIG. 6: (Color online) Variation of As-Fe-As bond angles with 
x at 1.7 K and with temperature for different K substitutions. 
The top panel shows ai, a' 2 and a'2 in orthorhombic setting. 
Blue circles represent ai, red triangles and black squares rep- 
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guides to the eye. 
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FIG. 7: (Color online) Neutron diffraction at 1.7 K with 
the magnetic Bragg peaks at a d-spacing of 2.45 Aand 
3.43 Aindicated by the arrows. They are absent above TV, 
for x = 0, 0.1, 0.21, 0.22 and 0.25. At x = 0.3, no magnetic 
peaks are observed. 

to the data. Removal of the time reversal symmetry from 
two of the mirror planes resulted in an antiferromagnetic 
arrangement of the magnetic moments with a magnetic 
wave vector Q = (1,0,1); that is, the Fe magnetic mo- 
ments are antiferromagnetically coupled in the x and z 
directions and ferromagnetically coupled along the y axis. 
This model is consistent with similar results previously 
reported for the parent BaFe 2 As2 material [15,16,44]. 

Rietveld refinements of both the atomic and magnetic 
structures were performed simultaneously as a function 
of temperature and doping, allowing the magnetic mo- 
ment to be defined in absolute units by normalization of 
the intensity of the magnetic Bragg peaks to the struc- 
tural Bragg peaks. The neutron data displayed in Fig. 
9 were collected at 1.7 K and normalized to the sample 
mass and exposure time (measured in beam pulses) . The 
figure qualitatively shows the intensities of the magnetic 
(121) and (103) Bragg reflections to remain roughly un- 
changed for the x = and 0.1 samples followed by a 
monotonic decrease upon increasing the K content until 
they nearly vanish at x = 0.24. The refinements show 
that the magnetic moment drops from \i = 0.75/ib for 
the parent BaFe2As2 material to 0.46/xb for x ~ 0.25 
(see Table 2). No magnetic peaks are observed beyond 
this value. 

The doping dependence of the two order parameters 
making up the AF/O phase is shown in Fig. 10. We 
compare 8 and /i 2 vs x, showing that they are directly 
proportional over the entire range of AF/O order. We 
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FIG. 8: (Color online) SQUID magnetization measurements 
for x = 0, 0.125, 0.15, 0.175, 0.2 and 0.24 in 2kG applied 
magnetic field. Insets are the first derivatives of the magne- 
tization curves, dM/dT (10 -5 ), used to determine the Neel 
temperatures given by the arrows. For x > 0.2, the magneti- 
zation anomaly at Tjv is too weak to be detected. 

have not been able to measure any samples between 
0.25 < x < 0.28, but a power law fit to 6 and /it 2 close 
to x c , i.e., to (1 — x/x c ) 2/3 , gives a critical concentra- 
tion of 0.252 with an exponent of /? = 0.125(1) for the 
structural and magnetic order parameters. The inset to 
Fig. 10 shows that TV is also proportional to fi 2 . In a 
mean- field model, Tjv scales as J^ 2 , where J is the effec- 
tive interionic exchange interaction, so this result would 
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FIG. 9: (Color online) Dependence of the neutron diffraction 
intensity for Bai_ a; K ;!; Fe2As2 at 1.7 K with x. The magnetic 
Bragg peaks are shown by the arrows. The solid lines repre- 
sent the calculated intensity of the Rietveld refinement. 
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FIG. 10: (Color online) Dependence of the square of the mag- 
netic moment, fx 2 (red stars), and the orthorhombic order pa- 
rameter, S = (a — b)/(a + b) (blue circles), at 1.7K on the 
potassium concentration, x. The inset shows a comparison 
of TV and ll 2 vs x, showing that Tjv oc fi 2 . Solid lines in 
the main panel and the inset are a fit to (1 — x/x c ) 2lS with 
/3 = 0.125. 

seem to indicate that J is approximately independent of 
x over this range. 



C. Superconductivity and Phase Coexistence 

Bulk superconductivity is observed for all samples with 
x > 0.15. Zero-field magnetization data shows that su- 
perconducting transition temperatures peaks at ~ 38 K 
for x = 0.4 before it slowly decreases to 3K for the end 
member KFe2As2 (Fig. 11). By comparing the samples' 
magnetic moment with that of a Sn-powder sample of 
similar volume, we estimate that these samples are bulk 
superconductors with a volume fraction of at least 80%. 
The uncertainty is due to variations in the demagneti- 
zation factor between samples. All the superconducting 
transitions are well-defined and sharp, even close to the 
critical concentration where T c is varying rapidly with 
x. Even low levels of compositional inhomogeneity asso- 
ciated with the uneven distribution of Ba/K ions would 
be revealed in magnetic susceptibility measurements by 
a broad or stepped-like transition from the normal state 
to the superconducting state, so this is further evidence 
of the sample quality. The increase in T c at low doping 
is approximately linear with dopant concentration up to 
x ~ 0.25, so we have estimated the critical concentration 
for the onset of superconductivity using linear regression 
to be x c = 0.130(3). 

In Rcf. [6] , we discussed the behavior of the order pa- 
rameters below T c . We observed a small reduction in 
both the magnetic and structural order parameters of 
approximately 5% at x = 0.21 and 0.24, without see- 
ing evidence of additional phases. In a scenario in which 
the sample divides into separate mesoscopic regions of 
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FIG. 11: (Color online) SQUID magnetization measurements 
(zero field-cooled) in 0.1 G magnetic field for (a) x = 0.15 
(solid squares), 0.175 (open squares), 0.21 (solid circles), 0.22 
(open circles), 0.25 (solid triangles), x = 0.3 (open trian- 
gles), and (b) 0.5 (solid triangles), 0.7 (solid circles), 0.9 
(solid squares) showing well-defined superconducting transi- 
tions. Magnetization values are normalized to mass of the 
samples, (c) Superconducting transition temperatures (onset 
T c ) of the underdoped compounds. Solid line represents the 
linear regression showing that the critical concentration for 
superconductivity is 0.130(3). 



AF/O phase and superconducting phase, this would im- 
ply that 95% of the sample remains in the AF/O phase 
and only 5% becomes superconducting. This is inconsis- 
tent with the magnetization measurements showing bulk 
superonductivity. While it is not possible to rule out the 
presence of other phases, the results indicate that there 
is microscopic phase coexistence of magnetism and su- 
perconductivity, with the reduction in the AF/O order 
parameters being due to a competition with the super- 
conducting order parameter. This competition has been 
discussed extensively by Fernandes et al [45], who show 
that there should be an additional phase boundary, with 
a positive slope vs x, between the coexistence region and 
the region of purely superconducting phase. We have not 
yet identified any anomalies corresponding to the phase 
line below T c , so we assume that it rises steeply with x. 
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FIG. 12: (Color online) Phase diagram of Bai_ a; K I Fe2As2 
with the superconducting critical temperatures, T c (circles), 
the Neel temperatures, Tn (stars), and the structural transi- 
tion temperatures, T a (squares). 

IV. DISCUSSION 

The overall phase diagram of Bai_ x K x Fe2As2 is shown 
in Fig. 12. We first discuss the nature of the spin-density 
wave order and orthorhombic order. Unlike the electron- 
doped compounds, where the two transitions split within 
increased doping, the two transitions are coincident and 
first-order in Bai_ a; K, r Fe2As2. We reported the first- 
order character of the transition by the observation of 
volume anomalies at T s [6]. Similar volume anomalies 
were also observed by Tegel et al [46] in unsubstituted 
SrFe2As2 and EuFe2As2 but, because of the small mag- 
nitude of these anomalies, the authors suggested that 
the structural phase transition may be second-order. 
However, other authors reported first-order transitions 
in polycrystalline SrFe2As2 [20] and single crystals of 
CaFe 2 As 2 [21] and BaFe 2 As 2 [44]. In the latter reference, 
a first-order-like hysteresis was obtained for the intensi- 
ties of the (101) Bragg peak when measured on cooling 
and warming. However, no such hysteresis was observed 
by Wilson et al [27], when examining their BaFe2As2 
single crystal. The systematic observation of volume 
anomalies across the phase diagram is unambiguous evi- 
dence that, at least in this system, all the transitions are 
first-order, although weakly first-order with extremely 
small hysteresis. 

Phenomenological theory of magnetoelastic coupling 
predicts the possibility of simultaneous first-order tran- 
sitions that are driven by a linear-quadratic term in the 
Ginzburg-Landau expansion [22,23] (see also Ref. [47]). 
This is the lowest-order term allowed by symmetry. If 
the magnetic transition were to occur at higher tempera- 
ture than the structural transition (in the absence of any 
competition), the magnetic order would drive the struc- 
tural order in a simultaneous first-order transition. The 



converse would produce two split transitions as seen in 
most of the iron-based compounds [24] . There is a report 
of a split transition in Bai-xK^F^Asa based on nuclear 
magnetic resonance results [25], but this is a local probe, 
which cannot necessarily identify compositional fluctua- 
tions. As we have already discussed, the neutron mea- 
surements, which represent true averages over the bulk, 
are quite unambiguous that the two transitions are si- 
multaneous, although there could be some rounding of 
the transitions at higher doping from small compositional 
fluctuations. 

The two order parameters are directly proportional 
to each other as a function of temperature, which 
seems to indicate an unusual biquadratic coupling in 
the Ginzburg-Landau expansion, rather than a linear- 
quadratic coupling. This is usually only observed at a 
tetracritical point [26,27] where two phase boundaries in- 
tersect, whereas our observations extend over a range of 
compositions. There are a number of possible reasons 
for this. The most intriguing and exotic idea is that the 
AF/O order parameters are both secondary to another 
order parameter and directly driven by it. This would be 
the case in, for example, valley density wave theory, in 
which a mother density-wave drives both the magnetic 
and charge-density- wave orders [48]. A second explana- 
tion is provided by the recent theoretical work of A. Nev- 
idomskyy [49] , which uses a microscopic Kugel-Khomskii 
model to produce a biquadratic spin-orbital term in the 
free energy. A subtle, but ultimately more conventional 
explanation, is that the coupling is linear-quadratic after 
all, but the proximity to a first-order transition produces 
a temperature dependence that is approximately equiv- 
alent to biquadratic coupling to first order [29]. This is 
in the context of a theory in which Ising-nematic order, 
produced by an itinerant model of Fermi surface nesting, 
drives the structural transition. Support for this expla- 
nation is provided by Fig. 10, where the doping depen- 
dence of the magnetic and structural order parameters 
at low temperature indicates a linear-quadratic coupling. 
Whatever the eventual explanation, it is clear that this 
result is key to understanding the nature of the normal 
state and the role of nematic order in the eventual super- 
conductivity. 

The strong coupling between AF/O order param- 
eters persists into the regime of phase coexistence 
with superconductivity. We have already argued that 
Bai_ x K x Fe2As2 is characterized by microscopic phase 
coexistence because mesoscopic phase separation would 
result in a significant decrease in the volume fraction 
of the AF/O phase below T c , The consensus in fa- 
vor of microscopic phase coexistence has existed for 
some time in the electron-doped superconductors, such 
as BaFe2_ a; Co a ;As2 [24], where the phase boundary be- 
low T c to a non-magnetic, purely superconducting region 
has also been identified. Theoretically, this behavior is 
consistent with unconventional s± pairing of the Cooper 
pairs suggesting that itinerant long range magnetism and 
superconductivity may coexist and compete for the same 
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FIG. 13: (Color online) (Top panel) Magnetic and structural 
phase diagram of electron-doped Ba(Fei_ ;c Co :!; )2As2 and hole- 
doped Bai_ a! K a ;Fe2As2 with the superconducting critical tem- 
peratures, T c (squares), Neel temperatures, TV (stars) and 
structural transition temperatures, T 3 (circles). The a;-axis 
is normalized to the charge carrier per iron atom. Data for 
the electron-doped side where the transition temperatures are 
represented with open symbols are taken from Ref [50]. The 
error bars for TV and T a values in the hole-doped side are 
within the symbols. The dashed line enveloping the super- 
conducting dome represents the Lindhard function taken from 
Ref [33]. (Bottom panel) Charge carrier dependence of the As- 
Fe-As bond angles for both electron- and hole-doping. Solid 
triangles represent the results of our neutron diffraction study 
at 1.7K for the hole-doped Bai_ x K :E Fe2As2. At this temper- 
ature one of the As-Fe-As angles splits due to orthorhombic 
distortion below x = 0.3. Therefore, we took the average of 
these two splitting angles. The As-Fe-As bond angle data for 
the electron doped side is taken from Ref [51]. Solid lines are 
guide to the eye. 

electrons [45] . However, the idea of microscopic phase co- 
existence was more controversial in Bai_ x K. I ;Fe2As2 be- 
cause of local probe measurements that seemed to indi- 
cate a phase separation into mesoscopic regions of mag- 
netism and superconductivity [30,31]. Since the most re- 
cent /^SR data are also consistent with microscopic phase 
coexistence [32], it appears that the earlier reports may 
have been due to compositional fluctuations close to the 
phase boundaries and that microscopic phase coexistence 
has now been confirmed. 

Finally, we discuss the electron-hole asymmetry in the 
phase diagram, shown in Fig. 13, where we have added 
data from the literature [50,51] to allow a comparison 
with the more commonly studied electron-doped super- 
conductors. In this phase diagram, the £-axis is normal- 
ized to the number of charge carriers per Fe atom. Neu- 
pane et al have recently suggested that this asymmetry is 
due to differences in the effective masses of the hole and 



electron pockets [33]. This is justified by ARPES data 
that show that hole doping can be well described within a 
rigid band approximation [52] . An ab initio calculation of 
the Lindhard function of the non-interacting susceptibil- 
ity at the Fermi surface nesting wavevector shows exactly 
this asymmetry, with a peak at x ~ 0.4 where the max- 
imum T c occurs. Our recent inelastic neutron scattering 
measurements of the resonant spin excitations that are 
also sensitive to Fermi surface nesting have shown a simi- 
lar correlation between the strength of superconductivity 
and the mismatch in the hole and electron Fermi surface 
volumes [34] , that is responsible for the fall of the Lind- 
hard function at high x. An overall envelope may be 
drawn (dashed line in Fig. 13) to encompass both the 
hole and electron superconducting domes of the phase 
diagram. If anything, the Lindhard function underesti- 
mates the asymmetry, predicting a larger superconduct- 
ing dome on the electron-doped side. We attribute this 
behavior to the fact that the iron arsenide layers remain 
intact in the potassium substituted series, whereas Co 
substitution for Fe disturbs the contiguity of the FeAs4 
tetrahedra and interferes with superconductivity in these 
layers. 

Interestingly, the maximum overall T c also correlates 
with the perfect tetrahedral angle of ~ 109.5° as demon- 
strated in the bottom panel of Fig. 13. In the plot, aver- 
age <As-Fe-As> bond angles for our K-substituted series 
have been extracted from the Rietvcld refinements. The 
As-Fe-As bond angles for BaFe2- a; Co. r As2 are extracted 
from the literature [51]. The continuity of the bond an- 
gles across the electron-doped and hole-doped sides of the 
phase diagram is remarkable and the crossing of the two 
independent angles at x ~ 0.4 to yield a perfect tetrahe- 
dron and maximum T c is clear. This has been remarked 
before in other systems [35,53]. It is possible that these 
two apparently distinct explanations for the maximum 
T c are two sides of the same coin. In a theoretical anal- 
ysis of the 1111 compounds [38], it has been suggested 
that the pnictogen height is important in controlling the 
energies of different orbital contributions to the cZ-bands 
and so affect the strength of the interband scattering that 
produces superconductivity. 

We now turn our attention to the SDW region of the 
phase diagram. While it is clear that spin-density-wave 
order has to be suppressed in order to allow supercon- 
ductivity to develop, it is not immediately clear what 
is responsible for the suppression. Both the strength of 
magnetic interactions and superconductivity, at least in 
an itinerant model, depend on the same Lindhard func- 
tion [54], the former on the peak in the susceptibility at 
the magnetic wavevector, and the latter on an integral 
over the Fermi surfaces. It would seem therefore that 
the magnetic transition temperature should also peak 
at x ~ 0.4. One intriguing reason why it would peak 
at x — is because magnetic order is more sensitive 
to disorder-induced suppression of the peak susceptibil- 
ity whereas superconductivity is more robust. There is 
some support for this idea from the observation that iso- 
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electronic doping produces a similar suppression of mag- 
netic order than seen with hole-doping [51]. On the other 
hand, Kimber et al succeeded in rendering the parent 
BaFe2As2 material to exhibit zero resistance at 30.5 K 
by the application of significant external pressures up to 
5.5 GPa [55], i.e., without introducing disorder, but they 
remarked that superconductivity needs to be confirmed 
by other bulk measurement techniques. Interestingly, the 
authors also correlate the induced T c with approaching a 
perfect tetrahedron angle of 109. 5o similar to our obser- 
vations for x <~ 0.4. 

In summary, we have synthesized high quality sam- 
ples covering the full phase diagram of Bai_ a; K a ;Fe2As2. 
Using high resolution neutron powder diffraction and 
SQUID magnetization measurements, we have investi- 
gated the effects of potassium substitution on supercon- 
ductivity, structural transformation and magnetic order- 



ing. Our measurements allowed the construction of a 
detailed magnetic and structural phase diagram, which 
displays a narrower phase coexistence region than the 
previous reports. Moreover, neutron diffraction and 
the SQUID magnetization data confirmed that magnetic 
and structural transitions are coincident with first or- 
der transitions. Additionally, we determined the effects 
of temperature and substitution on the various internal 
atomic and structural parameters. Our results confirm 
the importance of obtaining precise structural parame- 
ters across the whole phase diagram as a way of provid- 
ing insight into the nature of the phase competition that 
underlies iron-based superconductivity. 
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